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5h ! Abstract 

Fourth order derivative gravity in 3+1-dimensions is perturbatively renormalizable and is shown 
to describe a unitary theory of gravitons in a limited coupling parameter space. The running 
gravitational constant which includes graviton contribution is computed. Generically, gravita- 
tional Newton's constant vanishes at short distances in this perturbatively renormalizable and 
unitary theory. 
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1. Introduction and Theory 

(N 

J> " Newton's gravitational constant G in classical Einstein-Hilbert (EH) action is one of the 

smallest fundamental coupling constant with immense success in Cosmology |l|, |2|. Yet it is 
less understood quantitatively when we quantize the theory. The primary reason being that 
in 3 + 1 space-time dimensions this coupling constant behaves as (Mass)~ 2 and the theory 
has ultra-violet (UV) divergences. It was known Q that even quantum matter fields in the 
presence of background gravity demands that up to fourth order metric derivative terms need to 
be considered for UV regularization. Renormalization group (RG) studies of EH action [4( also 
suggest that we augment the gravitational action with fourth order derivative terms and witness 
UV renormalizability. Most general action functional up to fourth order derivative of metric 
is given by S, 
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where R is the Ricci scalar and R^ v is the Ricci tensor. The last term is also proportional to 
square of conformal Weyl tensor in 3 + 1 space-time dimensions, u) is dimensionless and M has 
dimension of mass. It was first shown in [J] that the action S is UV renormalizable using the 
4 — e dimensional regularization scheme. The primary observation being the action is non-linear 
in metric and hence has many number of basic interactions between gravitons. However, in the 
UV regime the fluctuating metric field (in 3 + 1 space-time dimension) has zero mass dimension 
and all these terms due to general co-ordinate invariance can be tamed by renormalizing the 
coupling parameters in Eq. (fTJ). 
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In the action S, classically if M 2 — > oo and A — > 0, only the EH term dominates, thereby 
accounting for all the known success of EH action. A deeper look into quantizing the action 
using ADM formalism showed that indeed S corresponds to the well-known constraints, 

space-time translation and a host of other constraints. The Hamiltonian itself reduces to a surface 
integral just as in pure EH gravity. Furthermore, the Hamiltonian is shown to be positive if the 
last term in S is absent and the coefficient of R 2 term is positive 0, Q . It is also shown that in 
Minkowskian background metric, the theory describes massless spin-2 gravitons, a massive scalar 
and a massive spin-2 field. 

Perturbation theory about a flat background <? M „ = rj^ + h plJ (?? M „ = {1,-1,-1,-1}) was 
studied in 0, [l(|. Here we investigate the action S for A = in the Landau gauge d^h pp = 0. 
The Feynman propagator in momentum space is, 
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where q is the momentum of fluctuating field /i M „. Various spin projectors are [Pi)^,^ = 
\ [T^oTup +T M pT va ] - \T pv T a p, (P s )^, a p = \T pu T a p, where T pv = r)^ - q p q„/q 2 . The first 
term of Eq. @ is the massless spin-2 graviton with two degrees of freedom, the second term 
is the scalar (Riccion) of mass Mj y/uj and the last term is that of massive spin-2 with mass M 
(M-mode), which is in agreement with the general ADM analysis. Furthermore, the residues at 
the pole for graviton and Riccion are positive, while for the M-mode it is negative. The M-mode 
renders the theory non-unitary. 

We now investigate the issue of the unitarity of the higher derivative gravitational action 
Eq. (TTJ) in the context of background Minkowskian perturbation theory. The main observation 
being that in quantum theory, all the couplings cease to be constant and will start running with 
momentum /i. Thus effectively at energy scale /x we have running couplings G(/i), M 2 (p) and 
This is an efficient way to incorporate the quantum effects in the renormalizable theory. 
In such a context if M 2 (/x)//i 2 > 1 for all /x 2 , then in any scattering process involving gravitons 
and Riccions we will never encounter the M-mode particle as an intermediate state to any order 
in Feynman loop expansion. Consequently the negative norm does not violate unitarity of the 
scattering matrix in graviton and Riccion subspace (physical GR). To reiterate consider a very 
high energy graviton scattering off with another graviton or Riccion. In principal there is enough 
energy to produce an M-mode, however the couplings and the masses are evolving with energy. 
In the renormalized theory where M 2 (/x 2 )//x 2 > 1, the mass of M- mode is always greater than 
the available energy, hence it is never realized on shell. In a scattering process of high energy 
graviton and Riccion going to other gravitons and Riccions with different momenta, M-mode 
can give off mass shell contributions, which affects the real part of the transition amplitude and 



never to the imaginary part (Cutkosky cut) Since M (p, > 1 implies it cannot go 



on-shell. Indeed we can reinterpret M-mode as an effective ghost normalizing the functional 
integral as conjectured in We will show that this can be explicitly realized in a suitably 
defined parameter space of the action S. 

2. Renormalization Group Analysis 

We first make the observation that our choice of parameterization in Eq. ([T]) naturally shows 
that the Feynman loop expansion is the same as perturbation theory in small G, namely G 
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effectively plays the same role as h. To one-loop the beta function of couplings have been 
computed using Schwinger-Dewitt technique 12fl . Using dimensional regularization in 4 — e 
dimensions and minimal subtraction [13l Il4j in the background Landau gauge, the RG flows are 
the following 15j: 
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where t = ln(/i//io) and all rhs contain the leading contribution in G (M 2 G is also taken to be 
small), with higher powers coming from higher loops being neglected. Using Eqs. ([3] and 2]) we 
solve for the running of ui, 
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where u)\ = 0.0457 and ui^ = 10.9343. Eq. ([7]) shows that — ui\ and —loi are two fixed points, 
with the former being repulsive and latter attractive under UV evolution or increasing t. As we 
see from our propagator Eq. ([2]) only u> positive can be in the physical domain. That is both the 
fixed points are in the unphysical domain. The rhs of Eq. ([7]) is always positive for u > 0. This 
means that a; is a monotonic increasing function of t and vice versa. Eq. readily allows us to 
express M 2 G in terms of t with which we can integrate Eq. (J7J to obtain, 
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where T = 133 ]{^ G and a — 399/50(^2 — uji), with subscript meaning that the coupling 
parameters are evaluated at t — or /i = fio- 

Eq. |J7J| transforms any evolution in t to evolution in u). Using Eqs. ([S]and[7]) we get, 
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G gets maximized at oj = -y/7/40, which lies in the physical domain. For convenience we 
choose this point to be our reference point /iq or t = and integrate Eq. (|9]) to obtain, 
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where Al = -0.3473 and A2 = -1.0027. Eqs. PandflO]! show that for large u or t, G ~ 1/w, 
thereby vanishing for large t 23], while for small uj, G ~ cj 7 / 20 reaching a peak at wo = y/T/40. 
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Similarly, using Eqs. (Q] and [5]) we obtain the running of M 2 /ia, which along with with Eq. ([7]) 
is integrated to give, 
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where Bl = 1.0802 and B2 = 0.2698. We can now analyze M 2 //x 2 from Eq. (JTT]). Alternatively, 
it is instructive to note that from Eqs. ([3j [5] and [7]) we obtain the running of ln(M 2 //u 2 ), 



ln / - j = vn_^o — ^ — hegz , (12) 
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This shows that M 2 //i 2 reaches a minima for w = oj* given by |u;» + ^ — ^g^j- 
Hence by demanding 
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we make the M-mode not realizable in the physical GR sector of the theory. The inequality is 
easily achieved by choosing // 2 G* appropriately. Perturbative loop expansion requires that M 2 G 
is small. Therefore M is a sub-Planckian mass, yet the running mass as dictated by interactions 
makes it physically not realizable even in post Planckian regime. 

Analogous to Eq. (TTT|) from Eqs. ((H [5] and [7]) we obtain, 

d - fM 2 \ 3 + ^ + ^, 
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showing that the Riccion mass relative to y, decreases monotonically. By a suitable choice we 
can make the Riccion to be physically realizable or not. So we conclude that there exists unitary 
physical subspace only with the gravitons or along with Riccions. 

We find that the solution of A in the RG flows is dictated by M 2 and G. From Eq. ((6]) because 
of M 4 term A increase in the UV region. These terms can be explicitly eliminated by adding 
fermion fields to our action. For example adding two spin-i Dirac fields with mass (5/4) 1 / 4 Af 
and M/y/2cu [16| (these additional fermionic fields can be interpreted as ghosts normalizing the 
functional integral). Consequently the other coefficients in Eqs. ([31 0] and [5]) do change in rhs 
making luq, u>i, u; 2 and shift slightly but all conclusions remains unaltered. Furthermore, in 
this modified theory the cosmological constant can be made to vanish for large t. 

We now look at the generality of the results. Adding matter fields to our action S modifies the 
RG flow Eqs. (OIU[5]and|6]). ^ n particular in rhs we will have terms which are dependent on mass 
m and couplings defined within the matter theory but independent of M and uj. For example to 



one-loop in matter fields these have been computed [16J, [17|, [15[ . Assuming m/M < 1, where m is 
a typical matter field mass, we notice that lu$, 0)%, w 2 , and w» all shift by a small amount and all 
our conclusion remain unaltered. We remark that in the absence of higher derivative terms i.e. 
the R 2 and conformal term, running of G is dominated only by matter fields and G increases in 
the UV regime [17| . In the renormalizable theory given by S, the Riccion and M-mode reverse 
this tendency and make it vanish. 

All our calculations are done in the Landau gauge. These calculations can in principal be 
done in other gauges, which have additional gauge parameters 14, 15| . They do not affect Eqs. 
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Figure 1: Running Newton's gravitational coupling. 

(OH]) and in Eq. ([5]) there are additional terms in rhs proportional to M 2 and M 2 /u>, but not 
M 2 uj. Consequently w\ and 0J2 are unaltered but uj shifts. It should also be noted that in all 
these gauge parameter choices, for large uj, G still vanishes as 1/w, while for small uj, G does 
vanish with a different power (as opposed to lu 7 / 20 in the Landau gauge). 

The mass of the Riccion M/y/uj is seen to be gauge parameter dependent, however in all these 
cases it has positive norm. We also looked at the unitarity issue in Prentki gauge (dihi^ = 0, 
where i — 1, 2, 3) [18|. In this gauge the Faddeev- Popov ghost is not propagating and we have 
exactly gravitons, Riccion with positive norm and mass M / yfoj and spin-2 Af-mode with negative 
norm and mass M, just as in Landau gauge. In addition there is a vector particle with positive 
norm and mass M (same as M-mode, RG flow equations have not been computed in this gauge). 
We see that in all these gauges the physical GR sector scattering matrix is unitary. 

Our RG analysis is at G = A = 0. However it is should be noted that higher derivative 
gravity action Eq. (fTjl was also considered at a non-trivial fixed point where both G and A are 
non- vanishing [3, Ho] . At this fixed point unitarity issue was considered [13, El at one- loop 
in a generalized de-Donder gauge using non-minimal regularization scheme. By doing Wilson's 
RG analysis on Euclidean functional integral, it is shown that the fixed point enjoys necessary 
spectral positivity properties over a finite range of parameter u>. 

We conclude that for large u, G ~ 1/w is a gauge invariant result. In the UV regime \/M 2 G 
is the coefficient of conformal term in the action and this is finite while R and R 2 coefficients 
are becoming large as w. R 2 being always positive makes functional integral better behaved as 
well. In addition the renormalized action enforces R not to fluctuate away from zero. In infrared 
the EH term R naturally is dominant along with cosmological constant. In infrared if A = (by 
tuning) the vanishing of G again enforces that R fluctuate around zero. 




3. Discussion 

The physical domain < to < 00 translates to a finite range in t, namely 

train _ ^ / ^2 ^0 + ^1 \ < J_ ^ ^ / ^0 + CJi \ _ t max 
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Figure 2: Af-mode ratio M/fi (top) and Riccion mass ratio M/fiy/uj (bottom). Dashed lines for M*//J,* = 25 
and solid lines for M t /fJ.* = 1.05. 

From one-loop analysis we see that there is no solution for to from Eq. if t > t max . For 
t < t?mn! w becomes negative i.e. M 2 /uj the Riccion mass square becomes negative signaling 
the instability of the vacuum. This infrared issue needs to be studied, namely the effects of the 
cosmological constant. To the leading order in Feynman loop expansion, at both extremes G 
vanishes as 1/u for t —¥ t max or w 7 / 20 for t — > t m in- 

Including inflationary scenarios the entire cosmological phenomenon ranges from as early as 
10 _42 s to the current age of the universe 10 18 s, a span of about 130 e-folds. In the higher 
derivative gravity theory we have essentially the following dimensionful parameters Go, Mo and 
[io (wo is set by the maximum of G). These can be traded for one dimensionful parameter Go 
and two dimensionless parameters t max and M*//i*. By choosing t max large enough for example 
about 100, we plotted G/G in Fig. Q] and M/(i and M/^/ujfj, in Fig. [2j For t max = 100, if 
M*//Xh» > 22.52 then both M-mode and Riccion are not in physical GR sector. In the Fig. [2] 
the dashed line plot represent such a case for M*//i* = 25. If 1 < M*//z* < 22.52, then Riccion 
is realized. In the Fig. [21 the solid lines represent an extreme case for M*//j,* = 1.05. From Fig. 
|5] we notice that Riccion is physically realizable in high energy scattering processes at most for 
about three e-folds before t max = 100. 

From Fig. [T]we note that the entire cosmology may span around the maximum of G accounting 
for sub-Planckian regime. Eventually in the post-Planckian regime > 10 19 GeV (10 _43 s), the 
gravitation coupling goes to zero at some finite energy ^ e * max (from Eq. ([5]) for uj > ^/7/40 
and M 2 G positive, G exponentially reaches zero and remains at zero). Naturally the question of 
physics beyond this scale is raised. 

To discuss the range of validity of this analysis, we rewrite the propagator Eq. ([2]) 
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From Eq. (p} and Eq. (IT6l) we see that there are three perturbative parameters G, M G and 
M 2 G and in the UV limit the important dimensionless perturbative parameters are latter two. 



G 



Under these circumstances we find M 2 G/oj and G both vanish at a finite energy scale /j,oe tmaa= . 
These results are valid if M 2 G is small as shown in the figures. However from Eq. (|3|) we also 
know that eventually M 2 G does become large before it reaches Landau singularity. We envisage 
two possible alternatives. Higher loop corrections may introduce terms like (M 2 G) lnw etc, such 
that the Landau singularity is avoided and t max approach infinity. Another alternative could be 
M 2 G becomes of order one, consequently we need to consider the conformal term in Eq. (TTJ) 
non-perturbatively. 

Finally we conclude that the action S, Eq. (TTJ) describes a perturbative quantum gravity as 
self consistent, renormalizable and unitary theory of gravitons and the curvature cannot become 
singular, in particular it cannot fluctuate wildly at sub-Planckian (22j or post Planckian regimes 
consistent with known cosmology. Its predictions asymptotically beyond Planck scale needs to 
be investigated further. 
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